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COMPARISON OF THEORETICAL AND EXPERIMENTAL ZERO-LIFT 
DRAG-RISE CHARACTERISTICS OF WING-BODY-TAIL 
COMBINATIONS NEAR THE SPEED OF SOUND 


By George Н. Holdaway 


SUMMARY 


The zero-lift drag rise at low supersonic speeds computed by 
linearized theory for several wing~body-tail combinations is compared 
with experimental data. The experimental data were obtained by the 
free-fall technique covering a Mach number range of М = 0.8 to M = 1.14. 
The procedures used in the theoretical calculations are outlined in 
detail. The results indicate that the computation method is capable of 
computing the drag rise at low supersonic speeds to an accuracy of about 
20 percent. 


Implications of the theory with respect to the selection of area 
distributions which will give a drag reduction over a range of Mach num- 
bers are examined. For cases where fuselage indentation is not feasible, 
* calculations indicate that drag reductions may still be possible by 
increasing the fuselage volume in front of and behind the wing. 


а 


INTRODUCTION 


Recent experimental results (refs. 1, 2, and 3) have demonstrated 
that the transonic zero-lift drag rise of wing-body combinations is 
primarily a function of the magnitude and rate of change of cross- 
sectional area along the longitudinal axis. This concept was utilized 
in the referenced tests to reduce the drag rise by indenting the body 
of a wing-body combination. These experimental results may be considered 
to be a qualitative verification of the linear theory as developed in 

н references 4 through 8. 


The purpose of this report is to examine the quantitative relation- 
> ship between the theory and experimental date from free-fall tests of 
several wing-body combinations. In addition, the implications of the 
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theory are examined with regard to modification of the area distribution 
of a configuration to keep the drag rise low over a selected range of 
transonic Mach numbers. 


This investigation utilized test data covering a Mach number range 
of M = 0.8 to М = l.ll} and а Reynolds number range of 2,500,000 to 
17,000,000, depending upon the wing mean aerodynamic chord of the con- 
figuration tested. 
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SYMBOLS 
dimensionless coefficients defining the magnitude of ihe 
harmonics of & Fourier sine series 


speed of sound, ft/sec 
drag at zero lift 


zero-lift drag coefficient 
g 2 абу 


zero-lift drag-rise coefficient, = 


zero-lift drag rise at transonic speeds, lb 
fuselage maximum diameter, in. 
fuselage or body length, in. 


Mach number, X 


number of terms or harmonics used in the Fourier sine series 
l 2 
dynamic pressure, 2 рү, 1b/sq ft 
projection of Sg on a plane perpendicular to x axis, sq in. 


cross-sectional areas formed by cutting the configurations 
with perpendicular or oblique planes X, sq in. 


wing area, sq ft 
free-stream velocity, ft/sec 


Gistance along the x axis measured from the midlength 
position, in. 


Cartesian coordinates as conventional body axes 
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0 angle between the z axis and the intersection of the cutting 
planes X with the yz plane 

Ц Mach angle, arc sin = 

p mass density of air, slugs/cu ft 

Фф transformation of the length x to radians, are cos иж 
radians 1/2 

X a series of parallel cutting planes intersecting the x axis 


at the Mach angle Н 
(At М = 1.0 these planes are perpendicular to the x axis; 
i.e., parallel to the yz plane.) 


у angle between the y axis and the intersection of the cutting 
planes X with the xy plane, arc tan [/M?-1 cos Ө 
g! first derivative of the projected cross-sectional area, = 
2 
5" second derivative of the projected cross-sectional area, 225 
THEORY 


Development of the Theory 


Contributions to the development of the theory were made by the 
investigations reported in references 4 through 8. The calculated drag 
rise of slender bodies of revolution and the analogy between the drag- 
rise equations and the induced-drag equations for a wing are presented 
in reference 1. The development of equations for bodies of revolution 
with minimum drag rise is presented in reference 5 which also illustrates 
the application of a sine series to evaluate the theoretical drag rise. 
The theory was extended to arbitrary systems in reference 6, which showed 
how these arbitrary systems could be represented by equivalent bodies of 
revolution. This latter theory was applied to lifting surfaces in 
references 7 and &. 


The application (and limitations) of these methods to the theoreti- 
cal computation of the drag rise of wing-body combinations is presented 
in reference 9. ‘The methods of reference 9 were used to make the theo- 
retical calculations for this report. The theory merely applies to the 
wave drag; any local shock or separation effects which might occur due 
to shape modification were not evaluated in the development of the theory. 
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Concepts Leading to the Drag Equation 


The derivation of the drag equation for a wing-body-tail combination 
is based on the theory that the configuration may be represented by a 
series of equivalent bodies of revolution. This theory is dependent on 
а simplified relationship between source strength and cross-sectional 
area. This relationship is used in planar wing and slender body of revo- 
lution problems. Specifically, &s was pointed out in reference 9, the 
source strength is &ssumed to be proportional to the normal component of 
the stream velocity at the body surface. The theory also assumes that 
the configuration is of a conventional type with thin symmetrical airfoil 
surfaces and a high fineness ratio body. Further exceptions and limi- 
tations to the theory are given in reference 9. 


The development of equivalent bodies of revolution will be illus- | 
trated by using the configuration shown in figure 1. The wing-body-tail 
combination is cut by a series of planes which always intersect the longi- 
tudinal axis at the Mach angie р. In other words, these planes are 
tangent to Mach cones. The plane identified in figure 1 as X; repre- 
sents one plane of a series of parallel planes which cut the configuration 
along the entire longitudinal axis. Each plane of this series intercepts 
the yz plane іп a line which forms the angle Ө, = 0° with the z 
axis. Similarly, planes Xə form the angles Ө with the 2 axis and 
the yz plane. For any one cutting plane of the series of planes 
Xo = f (05, н) the oblique cross-sectional area is projected on a plane 
perpendicular to the x axis. This projected cross-sectional area is 
plotted as a function of x. The resulting plot may be considered as 
representing the longitudinal distribution of cross-sectional area 8(х) 
of an equivalent body of revolution for the series of planes 
Xo = f (0, н). For any one value of u this process is repeated for 
other values of Ө ranging from Ө = О to 6 = 2л. However, if the con- 
figuration is symmetrical with respect to the xy and xz planes, then 
equivalent bodies for Ө from О to x only need be obtained. For 
bodies of revolution the area distribution is independent of Ө. 


With these concepts and with the use of the simplified relationship 
between source strength and cross-sectional area, the equation for the 
zero-lift drag rise &s & function of the rate of change of cross-sectional 
area can be derived from equation (46) of reference 7 and written as: 


pv? 21 nl i 
АР = — — n li S" (xi) S"( x5) in [xi = Xo |dxidxad6 (1) 
8х2 Jo o "o 


where ху and x» are two different locations along the x axis, 


S"(x) = — (2) 


Mo 


МАСА RM A53H17 „ТЕ, 5 


This generalized equation can be simplified by solving the double 
integral of the functions of x through a Fourier sine series in the 
same manner as used in reference 5. 


N 
ГА РА 2 
n n S"(x4) S"(x5) in] x, - Xp |dx1dx2 = -— z) ) nAn? (3) 
О О 


n=l 


(The factor 144 appears because the body length, 1, has been defined in 
inches for this report.) 


Where 
N 
ГА 
S'(x) DA H An sin n9 (4) 
п=1 
апа 
Ф = E 
arc cos 77 (5) 
then . 
x 
Аһ = 2 RI S'(x) sin пфаф (6) 
X 14, 


where the coefficients are a function of 0, since S'(x) is a function 
of 0. With this solution the simplified equation can now be written as 


l pv? ИМ” 
хо = р 8 (2) 8] Узы а i 


The computing procedure followed in applying the foregoing equations 
and theory to the determination of the zero-lift drag rise is presented 
in the Appendix of this report. 


CONFIGURATIONS AND TESTS 


Plan-view sketches of the models tested, and also the axial distri- 
bution of cross-sectional area normal to the longitudinal axis, are’ shown 
in figure 2. The different configurations will be referred to as models 
A, B, C, and D ав follows: 

Model A: aspect ratio 4 triangular wing with fuselage and tail 


Model B: aspect ratio 3 straight wing with fuselage and tail 


ents: 
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Model C: aspect ratio 6, 45° sweptback wing with fuselage and tail 


Model D: fuselage and tail (consisting of two vertical and two 
horizontal surfaces) 


General geometric data for all the models are presented in table I, with 
greater detail given for model D in figure 3. The fuselage and tail were 
the same for all models. The fuselage ordinates from the 8-inch to the 
139.4-inch station are given by the equation shown in figure 3. All 
cross sections of the fuselage were circular, and the nose of the fuse- 
lage was faired from the ordinate at the 8-inch station down to a 1-1/2- 
inch-diameter boom 15 inches forward of that station. The fuselage was 
not indented. 


The free-fall technique employed is described in detail in refer- 
ence 10. The tests covered a Mach number range of 0.8 to 1.14 and had 
approximately the following Reynoids number variation: 


Model M = 0.8 M = 1.1 


7,500,000 16,600,000 
5,750,000 12,400,000 
2,600,000 5,600,000 
35,800,000 77,000,000 


осо 


The values of Reynolds number for models А, В, and С are based on the 
wing mean aerodynamic chords, and the values for model D are based on 
the fuselage length of 210.5 inches. 


The accuracy of the measurements of the drag rise divided by q 

and the Mach number is believed to be within —2 = £0.06 and M = £0.01. 
Without an attempt to evaluate here the accuracy or limitations of line- 
arized theory, the theoretical computations as made are estimated to be 
accurate within 2 percent, except where linear theory might indicate a 
singularity in the drag curve. This accuracy of the theoretical data is 
based upon approximate methods for integrating the areas and determining 
the slopes of the area plots. Greater accuracies might be obtained by 
methods more exact than those outlined in the Appendix but probably are 
not justified. 


RESULTS 
The experimental zero-lift drag coefficients for the four models are 
presented in figure 4. The drag coefficients are based on the wing areas 


for models A, B, and C, and the fuselage maximum cross-sectional area for 
model D. The incremental drag rises above the subsonic values at M = 0.8 
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are presented in figure 5, together with the results of the theoretical 
computations. The rather large differences in wing area make comparison 
between configurations difficult; therefore, the data of figure 5 were 
multiplied by the respective areas used in computing the drag coef- 
ficients and the results presented in figure 6. The apparent late drag 
rise for model A is attributed to possible error in fairing the experi- 


mental data points which had scatter equal to 2 = 0.06 from the mean. 


DISCUSSION 


Agreement Between Experiment and Theory 


Comparison of the drag-rise values presented in figure 6 indicates 
that, in general, the theoretical computations give a good prediction of 
the experimental values. Limitations in the theory resulted in only 
qualitative agreement being obtained between theory and experiment at 
Mach numbers below about 1.02. At Mach number 1.0 the experimental drag 
rise is less than the theoretical value, varying from 46 to 74 percent 
of theory. It may be noted that the greater the departure of the con- 
figuration from that of a body of revolution the greater the difference 
between theory and experiment became for М = 1.0. 


For Mach numbers greater than 1.0, the effect of theory limitations 
is diminished and the level of the drag rise is predicted quantitatively. 
In general, the values above M = 1.02 are in error by about 20 percent 
with а maximum deviation of experiment from theory of 26 percent. 


Application of Theory 


The prior section illustrated the adequacy of the theory in estimat- 
ing the zero-lift drag rise of wing-body-tail combinations. It is of 
interest next to examine the implications of the theory with regard to 
possible drag reductions through modifications of the area distribution. 


Consider a design problem involving a configuration similar to 
model B with an engine or other components within the fuselage which 
might make fuselage indentation impractical. In a case such as this the 
drag reduction would have to be attempted by adding volume before and 
behind the wing. That such an approach might be successful is illus- 
trated in figure 7(a). The solid curve is the computed drag rise for the 
unmodified configuration. The lowest, straight line represents the theo- 
retical drag rise for a Sears-Haack body with the same maximum cross- 
sectional area as the wing-body combination. The large difference 
between these two curves suggests the possibility of reducing the drag 
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rise of the original configuration by modifying its area distribution 
to approach, as & limit, the Sears-Haack distribution. 


The modification curves presented in figure T(b) represent various 
fuselage revisions, obtained by adding volume to the fuselage in front 
of and behind the wing. Modification 1 is similar to the type of fuse- 
lage revision discussed in reference 1, in that it is based entirely on 
the distribution of cross-sectional area perpendicular to the x axis. 
Such & revision would be expected to reduce the drag rise &t & Mach 
number of 1.0, and this anticipation is realized in figure 7(а). However, 
the desired equivalent body of revolution for M - 1.0 is not the desired 
Shape for other Mach numbers, and & rspid rise in drag with increasing 
Mach number therefore results as shown in figure 7(a). This drag rise is 
independent of possible separation effects which were not evaluated in 
the theory used. 


Clarification of the above point is afforded by consideration of 
figure 8(a). On this figure are shown the projected cross-sectional 
area, distributions (S) of the original wing-fuselage-tail configuration 
for various angles \Ущах. The largest angle Ymax is used for each 
Mach number computation and is the complement of the Mach angle. (The 
equivalent tail area for the larger angles showed only slight differences 
from the Ymax = 8.089, во was omitted from this figure for clarification. 
The fuselage having a high fineness ratio had essentially the same area 
distribution for the entire Mach number range.) Shown on figure 8(a) is 
also а Sears-Haack body area distribution for a body with maximum cross- 
sectional area, equal to that of the original wing. plus body. Modifi- 
cation 1, which optimizes the area distribution on the original wing- 
fuselage combination for a Mach number of 1.0, consists of additions to 
the body cross-sectional area sufficient to fill in the difference indi- 
cated by the shaded area. The fact that excess area is added for other 
Mach numbers is indicated by the penetration of the (S) curves for 
other values of Vmax into this shaded area. If it is desired to reduce 
the drag at Mach numbers somewhat above 1.0, much less volume should be 
added to the regions immediately in front of and behind the wing. 


In arriving at such a compromise for the wing body under consider- 
ation, the following approach was followed to establish the lines for 
modification 2. The five area distributions (S) shown in figure 8(a) 
were arbitrarily replaced by & single curve representing the arithmetic 
average of the five curves. The Sears-Haack body area distribution was 
modified to represent a Sears-Haack body having & maximum cross-sectional 
area equal to the maximum of this averaged curve. The volume added to the 
fuselage was determined from the difference between this new Sears-Haack 
body curve and the average curve as shown in figure 8(b). The Ypax 
curves for М = 1.0 and М = 1.1 are also shown on figure 8(b) to illus- 
trate the reduction in the size of the new Sears~Haack body as well as 
the reduction in volume to be added immediately before and behind the 
wing. The computed drag rise for this second modification (fig. 7(a)) 
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shows an improvement in comparison to the original unmodified configur- 
ation over the entire design Mach number range of from 1.0 to 1.1}. 


Modifications 3 and 4 were derived by the method described in 
reference 9 for the design of a fuselage modification optimized for a 
specific Mach number. Mach numbers of 1.05 and 1.14 were arbitrarily 
selected and the results of the calculations are presented in figures 7(a) 
and T(b). Modification 3 produced only a slight improvement over modi- 
fication 2, but it was maintained over the Mach number range shown. 
Apparently, as the minimum drag curve for a configuration is approached, 
slight streamwise modifications to the fuselage about the optimum shape 
can be made without greatly changing the drag. Further reduction of 
the drag rise at the higher Mach numbers at the expense of an increase 
at the lower Mach numbers could probably be effected by optimizing for 
the highest operating Mach number as was done in modification №. Due 
to the limitations of the drag computations near Mach number 1.0; 
experimental investigations should be made to see if these small drag 
differences for modifications 3 and } really do exist. 


CONCLUSIONS 


The results of the computations and experimentation presented in 
this report indicate that for the models tested the following conclusions 
are justified: 


1. А computation method has been established which is capable of 
predicting the zero-lift drag rise at low supersonic speeds for a wide 
variety of wing-fuselage-tail combinations within an accuracy of about 
20 percent. 


2. In establishing the modifications to be made to a fuselage to 
reduce the drag rise, the effect of various revisions should be examined 
for & range of Mach numbers. 


3. For cases where fuselage indentation is not feasible, theory 
indicates that drag reductions may still be possible by increasing the 
fuselage volume in front of and behind the wing. 


Ames Aeronautical Laboratory | 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 17, 1953 
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APPENDIX 
PROCEDURE USED FOR COMPUTING ZERO-LIFT DRAG RISE 


The detailed application of the drag equation developed in the 
theoretical section of this report is outlined in the following para- 
graphs. The general procedures are presented with various simplifi- 
cations and suggested methods based on experience gained in making the 
computations for this report. The simplest case is presented first in 
order to better illustrate the basic steps and to shorten the presen- 
tation of the more involved procedures. | 


Calculation at М = 1.0 for both Symmetrical and 
Unsymmetrical Configurations 


At M= 1.0 the cutting planes X are perpendicular to the x 
axis and the projected cross-sectional areas are independent of the 
angle 9. Therefore, the drag rise at M = 1.0 is computed in the same 
manner for symmetrical or unsymmetrical configurations. | 


Step 1.- Determine the total cross-sectional area of the wing-body- 
tail combination and plot the area distribution as shown in figure 9(a). 


Step 2.- Differentiate the area plot, figure 9(а), with respect to 
x, and plot the results as a function of Ф as is indicated in 
figure 9(b). The technique of reference 11 may be used to differentiate 
the area curve. 


Step 3.- Fit a Fourier sine series to the curve of figure 9(b) in 
order to determine the magnitude of the coefficients of the various har- 
monies. For example, from equation (6) of the theoretical section the 
coefficients An are computed as follows: = 


x 
1 zi 25 sin 989 Е 
xi 
e 
» Pas : 


The integrations may be accomplished by pünch card machine computation, 
electronic-wave-analysis equipment or by ordinary long-hand computing 
methods. However, early investigations in connection with the compu- 
tations made for this report showed that a large number of terms (over 12) 
would be required to fit the irregular curves normally obtained, which 
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would rule out long-hand computing methods. Available electronic-wave- 
analysis equipment gave the relative magnitude of the coefficients but 
not their absolute values. The punch card computing methods were used 
for this report. Practical considerations of over-all accuracy did not 
justify calculating for more than 24 harmonics. On the other hand, com- 
puting for less than 12 harmonics reduced the resultant theoretical drags 
as much as 40 percent for some of the configurations investigated in 
this report. 


Step 4.- Compute the drag rise from equation (T) which for М = 1.0 
reduces to the following coefficient form: 


n-224 


2 
C ———) z 2 
Ао = 576 S. NP "^n 


Il 


In using equations (6) and (7) they may be simplified by deleting 2/1 
and (1/2)?, respectively, provided A, is considered ав a dimension in 
inches. 


Calculation for Mach Number Range, Configuration Symmetrical 
With Respect to xy and xz Planes, Wings 
and Tail Surfaces in xy Plane Only 


The first three steps in this calculation are essentially the same 
as for the М = 1.0 calculation except the computations should be made 
for at least four equivalent bodies of revolution rather than for just 
one. 


Step 1.- Select the Mach number range to be covered by the compu- 
tations and determine the cross-sectional areas of the equivalent bodies 
of revolution for each value of Ө (u = constant) and plot as shown in 
figure 9(a). By judicious selection of cutting planes for the compu- 
tation of the drag for the highest Mach number, the drag at lower Mach 
numbers can be computed from the same cross sections. This means that 
the values of Ymax used for the lower Mach numbers are used for the 
intermediate values of Ў for the higher Mach numbers. For thin air- 
foil surfaces, which for purposes of this computation may be assumed to 
lie in the xy plane, the surfaces may be cut at angles V by planes 
perpendicular to the xy plane. The actual obliqueness or error in the 
Wing cross sections disappears as the areas are projected onto the yz 
plane. It should be recalled that \ is the intercept angle of a Ө 
plane with the xy plane. Тһе prior statements can be illustrated by 
an example showing the angles Ymax which might be used in determining 
the cross sections for computations for M = 1.30. 


12 epum NACA RM A53H17 
Vay = arc tan JM -1 


M V max» deg 


1.00 О 

1.02 8.08 
1.02 11.38 
1.05 17.75 
1.15 28.70 
1.30 39.71 


The equivalent bodies of revolution are determined for each of the above 
angles by cutting the configuration by a family of parallel planes. Then 
the area for each cross section is projected onto the yz plane by multi- 
plying the area by cos шах. This produces the equivalent bodies of 
revolution used in the Pinal computations, where the drag computations 

for М = 1.05 will use only the equivalent bodies for the first four 


angles. 


For high fineness ratio fuselages, it is not necessary to cut the 
fuselage at various angles, because the projection of the cross-sectional 
areas onto the yz plane results in essentially a constant area distri- 
bution with Mach number as mentioned previously. For a fixed area distri- 
bution, the shape of the cross section of a high fineness ratio body is 
relatively unimportant if the body ends in a cylinder or a point (see 
ref. 12). The zero-lift drag-rise coefficients for the Sears-Haack 
bodies can be computed from equation (14), reference 5, which illustrates 
that the supersonic drag is constant: 


2 
EHO 


Step 2.- Differentiate the area plots and plot as a function of Ọ 
as indicated in figure 9(b). 


Step 3.- Fit a Fourier sine series to the curves of the area plots 
in the same manner as for the M = 1.0 computations. 


Step 4.- For the largest Mach number and angles Ў, compute the 
angles Ө and plot EnAy,“ as shown in figure 9(c). The same data can 
be used for the intermediate Mach numbers as follows: 
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“Angles Ө are taken symmetrical about Ө = 1/2. 


Step 5.- Integrate the area under curves similar to figure 9(c) and 
compute the drag-rise coefficient from the following equation derived 
from equation (7): 


ET ы 


nix Dm Su G 2) / a 


Calculation for Mach Number Range, Configuration Symmetrical 
With Respect to xy and xz Planes, Wings 
and Tail Surfaces in Both Planes 


The essential differences from the prior method is that the verti- 
cal surfaces are rotated 90° into the xy plane, and if the same values 
of Y used for the cutting of the horizontal surfaces are used to cut 
vertical surfaces then the cross sections for the vertical surfaces will 
correspond to different angles @ than those used for the horizontal 
surfaces. Because the areas must be combined to give an equivalent body 
of revolution for one value of Ө, cross plots of the areas should be 
made or the vertical surfaces ahouid be cut at different angles than the 
horizontal surfaces. A satisfactory and simple procedure is to cut the’ 
vertical surfaces for one or two angles more than the horizontal surfaces 
to ensure a uniform variation of angles from O to 909, and then proceed 
as follows: 


Step l.- Multiply the cross-sectional areas for the vertical 
surfaces by cosY. 


Step 2.- For each fuselage station plot the areas from step 1 as a 
function of Өү -.0g - 90°. Subscripts V and Н refer to the verti- 
cal and horizontal surfaces, respectively. 


«=== 
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Step 3.- Read the areas from step 2 that correspond to the values 
of Ө. 


Step 1.- Combine the cross-sectional areas (already multiplied by 
сов V for the vertical and horizontal surfaces to give one equivalent 
body for one angle Og = Ө. 


Step 5.- Compute the drag-rise coefficients as before for the case 
where all surfaces were in the xy plane. 


Calculation for Mach Number Range, Unsymmetrical Configuration 


This calculation is similar to the prior method except that the 
final integration must be performed from Ө = О to Ө = 2х and the 
equivalent bodies of revolution determined from negative angles of Y 
as well as from positive angles. These factors lengthen the problem 
but do not add to the complexity. 


Time to Perform Calculations 


Computations for this report, as was mentioned previously, were 
made by punch card machine computations. Time to compute the drag rise 
for a wing, similar to the models tested for this report, in combination 
with a high fineness ratio body from М = 1.0 to М = 1.2, was found to 
require approximately 80 computer hours. This value is based on experi- 
ence gained in making the computations for this report and includes the 
time to lay out the cross sections and integrate the areas graphically. 
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TABLE I.- DIMENSIONS OF TEST MODELS 


(a) Wings! 


Area, sq ft 
Mean aerodynamic 


7.36 
6 


Taper ratio 


Airfoil section 


Fineness ratio . . . . «e «© «© e © e © © © © © ss э э э ng] 
Maximum diameter, ft * * ^ е . е a е . е е LÀ * е e * . e Ф l . 417 
Nose boom diameter, in. Е > o . . 1.5 


Horizontal surfaces 


Area, ва ft 2. . 6 ee ee ee ee ee ee m om th n n t ng 6.0 
Aspect гаййо........... 2 om om ot m n t! t] n t|: 1.5 
Taper ratio . . ea e © © © © © © новое и 0.20 
Airfoil section ........... МАСА 65006 streamwise 
Sweep of 0.25 chord . . « eee eeso we 499. +5 4.50 


Vertical surfaces 


Area, ва ft оо 
Aspect ratio . . . « s « «© e 4 45 ж t t ot э э 91 n tg Dl 
Taper ratio . . s.. ee eo e eo a ae 0.22 
Airfoil section . . . . . NACA 650 to 0.25 chord 
Sweep of 0.25 chord А "M 450 


O9 perpendicular 


1See figure 2. 
“See figure 3. 
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Equivalent body of revolution 
formed by the planes Xo 


18459 


Figure 1.- Illustration of the cutting planes X and the angles 
Ф and V, which are the intercepts of these planes with the 
yz and xy planes, respectively. 
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SRE” 


Station 
0 
Plan view Station 
-25.1 
45° 
Angle-of-attack // | 
Vase ру rÍ = 62.5 
-«— — . A. —————— — 
Equation for body ordinates between x = 8 and 6 
= 139.4 igs: р = 8. i4 -1) 
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Note: All dimensions are in inches, “RA 


Figure 3.- Geometry and dimensions of model D. 
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Zero-lift drag coefficient, Cp, 


© 


Galen Tt NACA RM A53H17 


—— А 30.10 
UBRE 21.68 | Sw 
9.02 


————D 1.58 (Cross 
Sectional) 


| BEENENE 
„8 1. ‚2 


Mach number, М 


Figure }.- Experimental zero-lift drag coefficients, based upon wing 


area for models A, B, and € and upon maximum cross-sectional area 
for model D. 
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—— Experiment 


Zero~lift drag coefficient rise, ACDo 


O 
‚8 9 1.0 3.1 1.2 
Mach number, М 


(d) Model D. 


Figure 5.- Comparison of experimental drag rises with values calculated 
utilizing linearized theory. 
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(d) Model D. 


Figure 6.- Comparison of theoretical and experimental increases in drag 
divided by dynamic pressure. 
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Unmodified 
Modification 1 
Modification 2 
7 — Modification З 
— — — ~ Modification 4 
——— —- = Sears-Haack body 
with game S8max. 


ACpg 


Calculated zero-11ft drag rise, 
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Mach number, М 


(а) Calculated drag rise, Абу. 
Figure 7.- Comparison of the computed drag rise of model В with the computed values for the ваше 
configuration with fuselages of greater volume modified in various ways in an attempt to 
reduce the drag rise: 
(1) Enlarged body indented for equivalent wing area at М = 1.0 
(2) Enlarged body indented for the arithmetic average of equivalent wing areas (5) from 
М = 1.0 to M =1.14 
(3) Enlarged body optimized for M = 1.05 
(1) Enlarged body optimized for M = 1.14 
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Figure 7.- Concluded. 


Body stations, inches 
(b) Fuselage revisions. 
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Fuselage stations, inches 


(a) Modification 1 for M = 1.0. 
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to the fuselage ahead of and behind wing. 


Figure 8.- Modifications to the area distribution of ‘model B to reduce drag by adding volume 
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(b) Modification 2 for M = 1.0 to 1.14. 


Figure 8.- Concluded. 
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Figure 9.- The three steps in the theoretical solution of the drag rise. 


For an unsymmetrical configuration, plot C should be made from 


Ө = O to 8 = 2п 
ee es 
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